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1. Introduction 

Over many years, lattice QCD simulations have been done mainly for zero and two flavors of 
fermions because of the algorithmic difficulty of simulating odd flavors. Recently, though, several 
algorithms have been developed for odd-flavor simulations. A new comer is the rational hybrid 
Monte Carlo algorithm (RHMC) in which a fractional power of the fermion matrix needed to 
simulate odd fermions is given by a rational approximation [[T[]. The advantage of the RHMC is that 
the approximation error can be made small with a low approximation degree. 

One-flavor QCD had not been seriously studied. However, motivated by theoretically interest- 
ing properties ( one-flavor QCD has no pion and is expected to have no chiral symmetry ) a study 
of the one-flavor spectroscopy has now been started!^. 

One-flavor QCD at finite temperature was studied some time ago[^,^. Alexandrou et al. used 
the multi-boson algorithm to simulate one-flavor QCD and mapped a rough phase diagram of the 
one-flavor QCD in the heavy quark region. In order to locate the end-point accurately, they also 
used the effective 3D Potts model and estimated JQ. ~ 0.08[^. 

In this study we use the RHMC to simulate one-flavor QCD and try to locate the end-point by 
the Binder cumulant of the Polyakov loop norm. 

2. One-flavor algorithms 

The lattice QCD partition function is given by 



where D{U) is the fermion matrix, nf the number of flavors and Sg{U) the gauge action. We use 
the standard Wilson fermion and the standard Wilson gauge action. With this partition function the 
expectation value of an operator Q. is given by 



For multiples of even-flavor the fermionic determinant ( here for 2 flavors as an example ) is ex- 
pressed as 



the fermionic action, Sf = <p^{D^D)^^(p is manifestly positive. This positive form of the action can 
not be easily realized for odd-flavors. For this positive fermionic action one can easily update (p 
using the heat-bath method, i.e. <p = D^rj, where rj is drawn by P{ri) ~ e^^ ^ . 

If we use the identity of detD"' = UftrlogD we can make simulations of any numbers of 
flavors using the R-algorithm[^ which is, however, an inexact algorithm introducing errors to the 
results. When we calculate the fermionic forces in the R-algorithm, the calculations of tr{D^^D') 
appear. Such calculations are numerically very costly and usually they are estimated approximately 
using the random noise method, which introduces the approximation errors. These errors can not 
be removed completely. The order of the errors is 0{dt^), where dt is the step-size. In order to 
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(2.2) 
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obtain the exact results from the R-algorithm, one needs to extrapolate them to the zero step-size 
limit, i.e. 5t 0. 

The exact algorithms of odd flavors can be constructed using Liischer's idea[^. He approxi- 
mates the inverse of the fermion matrix using the following polynomial. 

i«n;Li(£>-zO, (2.4) 

where Zi are the roots of the polynomial. With this polynomial a multi-boson algorithm which 
is algorithmically very different from the HMC can be constructed. Originally it was applied for 
even flavors and later generalized to any numbers of flavors!^, ^]. Using the 75 hermiticity of D, 
D = 750^75, the determinant of a single D for «y = 1 is written as 

detD « Ati[Tli^{D)T„,2{D)]-\ (2.5) 
oc j Y^ll^<^D^te-^^^'^-'^^-'^^', (2.6) 
= Jw;/^D^jD^,e-^^^-^>, (2.7) 

where Tn/2{D) = n"^^(D — z,) and r\i = {D — Zi)<^i- Each bosonic action {D — Zi)"^ {D — Zi)(j)i is 
now positive and are updated using the heat-bath method. 

The same idea of using the polynomial can also be used in the framework of the HMC algo- 
rithm, de Forcrand and Takaishi first used Liischer's polynomial for the HMC in order to reduce 



the computational cost of the algorithm! 10]. With the help of the polynomial, an odd-flavor HMC 

algorithm is constructed [11] as follows. Using the 75 hermiticity of D, the determinant of D is 
written as 

detD « det[r+2(£>)r„/2(D)]-\ (2.8) 

oc |D0tD0e-^'[^"/^(^)^"/2(°)]^ (2.9) 

= Jo^'^D^e-'^''^, (2.10) 

T«/2/ 



where T„/2{D) = IT. ^^(D — z,) and rj = r„/2(D)0. Thus, as in the conventional HMC algorithm, 

n/2 



the fermionic action Sf = 0^[rJ,2('^)^«/2('O)]</' is made to be positive and is updated using the 



heat-bath method. 

Liischer's polynomial is not nessesarily the optimal one. Instead, one can find a rational ap- 
proximation with more accuracy. Clark and Kennedy constructed the HMC algorithm with a ratio- 
nal approximation [||], which will be described briefly in the next section. 

3. Rational Hybrid IMonte Carlo 

In the RHMC the fermion determinant is rewritten as 

detM« oc J D(l)'^D(l)e-^'^'"^ = J D0tD0e-^''•('^)"^ (3.1) 
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where a = nf /2, M = {D^D) and the rational approximation r„{M) ' is given by 

In this way, the fermionic action Sf = <p^rn{M)'^<p is made to be positive. Thus the framework of 
the HMC algorithm can be used with this approximation. 

The operations such as r(M)^'T] including (M — jS,)^'!] can be done by the multi-shift solver[ 12]. 
The computational cost of this operation is expected to be similar to the HMC [IT]]. 



The coefficients of the rational approximation are numerically calculated by the Remez method[]13p. 
The rational approximation error can be made small with decreasing order n. Typically, 0(20) or 
less gives the maximum relative error that is good to machine precision[|l|]. Fig.|] shows an example 
of the relative errors calculated on 20^ x 4 lattices. We calculate | r„ (M) " ^ — r,,,,,^^ (M) " ' p / 1 r,,,,,^^ (M) " ^ 1 
as a function of n. Here is random gaussian noise vectors and rimax is set to 60. The results shown 
on Fig.l are averages taken over 10 configurations at j3 = 5.63 and K = 0.12. The error bars are 
smaller than the symbols. We see that with increasing n the relative error decreases quickly. The 
relative error with « = 20 is already very small, less than 10^^^. For our simulations we take n = 20 
or 25( for the largest lattice 20^ x 4 ). 

Here, note that the low approximation degree itself does not mean that the cost of the RHMC 
is small. The rational approximation in the RHMC contains solver calculations such as D^^<p. 
Roughly speaking the cost of the RHMC is proportional to that of the solver calculations. On the 
other hand, in the polynomial HMC the calculation of is replaced with the calculation by the 
polynomial approximation. Therefore, the cost of the polynomial HMC is directly proportional to 
the number of the approximation degree. 

Fig.^ shows the plaquette values obtained from «/ = 1 simulations on a 6"* lattice at j8 = 
5.45 and K" = 0.160 as a function of the approximation degree. The results from the RHMC are 
consistent with those from the R-algorithm and the polynomial HMC unless the approximation 
degree is very small, like « ~ 5. 
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Figure 2: Plaquette values calculated on a 6^ lattice at j3 = 5.45 and K" = 0.16 as a function of the approxi- 
mation degree. The results from the polynomial HMC and the R-algorithm are taken from WM- 



4. One-flavor simulations 

We simulated one-flavor QCD on x 4 lattices with L = 8, 12, 16 and 20. We used an im- 
proved integrator, 2MN integrator, for molecular dynamics simulations in the RHMC. The 2MN 



integrator is 50% faster than the conventional 2nd order leap frog integrator [[14|]. The step-sizes 
are set to the values that give acceptances from 60% to 70%. These acceptances are shown to be 
optimal for any 2nd order integrator! 15]. We generated about 2 x 10^ (1 x 10^) trajectories for L=8 



(12~20). 

Fig.^ shows the Polyakov loop norm for each K" as a function of j8. At large K we see no 
evidence that on larger lattices the discontinuity across jSc will be pronounced. On the other hand, 
at small k the discontinuity appears for larger lattices, which may show evidence of first order 
phase transition. 

We calculate the Binder cumulant B4 defined by 

^4 = ^. (4.1) 

64 takes 3 or 1 for the first order phase transition or crossover. At the critical point B4 takes a certain 
value dependent on the universality class. For 3d Ising universality B4 is 1.604[|^]. Fig.^ shows S4 
for the Polyakov loop norm as a function of k. The results are still noisy and the clear critical point 
where all curves intersect is not determined yet. However we estimate it to be jq. = 0.07 ~ 0.08. 



5. Summary 

The one-flavor QCD simulations are done by the RHMC. The relative error of the rational 
approximation can be made small using a low approximation degree such as n ~ 20. We compare 
the plaquette results from the RHMC with those from the polynomial HMC and the R-algorithm. 
The results are consistent with those from the polynomial HMC and the R-algorithm. 
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Figure 4: Binder cumulant of the Polyakov loop norm for various lattice sizes as a function of K. 
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We calculate the Binder cumulant B4 of the Polyakov loop norm on x 4 lattices with L = 
8, 12, 16 and 20. The B4 curves are expected to intersect at the critical point Kc and we estimate 
roughly K-^ ~ 0.07 - 0.08. 

In future in order to locate the end-point accurately we plan to use a larger lattice and to 
improve our data statistics. 
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